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(N ; Abstract 

-(— > ■ 

O ■ Superfluid helium consists of two inter-penetrating fluids, a viscous normal fluid and an inviscid 

o 

superfluid, coupled by a mutual friction. We develop a two-fluid shell model to study superfluid 
turbulence. We investigate the energy spectra and the balance of fluxes between the two fluids 
1-^ ' as a function of temperature in continuously forced turbulence, and, in the absence of forcing, 

the decay of turbulence. We furthermore investigate deviations from the k~^'^ spectrum caused 
by the mutual friction force. We compare our results with experiments and existing calculations. 

c/j I We flnd that, at sufficiently low temperatures a build-up of energy develops at high wavenumbers 

-)-* ■ 

C^ ' suggesting the need for a further dissipative effect, such as the Kelvin wave cascade and plionon 

B 

I . emission. 

a, 

O . PACS numbers: 67.25.dk Vortices and turbulence in superfluid ''He 

47.32.C- Vortex dynamics 
^ ' 47.27.Gs Homogeneous isotropic turbulence 
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I. INTRODUCTION 



Helium II (the low temperature, quantum phase of liquid ^He) consists of two inter- 
penetrating fluid components [jj: the inviscid superfluid (associated with the quantum 
ground state) and the viscous normal fluid (consisting of thermal excitations). Each fluid 
component has its own density and velocity field, p**, u^ for the superfluid and p", u" for the 
normal fluid, where p = p* + p" is helium's density and the relative proportion of superfluid 
and normal fluid depends on the absolute temperature T. 

What makes helium II particularly interesting is that the superfluid vorticity is con- 
centrated in thin, discrete vortex filaments of fixed (quantized) circulation k, = h/m = 
9.97 X 10~^ cm^/s, where h is Planck's constant and m the mass of one ^He atom. The 
vortex filaments interact with the thermal excitations, causing a mutual friction between 
superfluid and normal fluid 2|. 

Turbulence in helium II (quantum turbulence) is easily generated in the laboratory by 
stirring or pushing liquid helium with grids and propellers |3|-l5| , vibrating forks 6| or other 
means. Quantum turbulence is a complex state [7| in which normal fluid eddies of arbitrary 
shapes and strengths interact with a tangle of discrete vortex filaments of fixed strength. 
Despite the two-fluids nature of quantum turbulence, experiments have shown remarkable 



similarities with turbulence in ordinary 
pressure drops alo ng p ipes and channels 



iuids (classical turbulence), for example the same 
8|, drag crisis behind a sphere J9|, and Kolmogorov 
energy spectrum [10|. The last property, which is our particular concern in this paper, 
describes how the energy of the flow is distributed over the length scales in homogeneous 
isotropic turbulence. 

Whereas turbulence in ordinary fluids is studied on the firm ground of the Navier- Stokes 
equation, there is not such an established equation for turbulent helium II. Existing numer- 
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ical simulations [Hi, Il2| of quantum turbulence follow the approach of Schwarz [13[, who 
calculated the properties of the vortex tangle by integrating in time the motion of a great 
number of individual vortex filaments in the presence of a prescribed normal fluid, neglecting 
the back reaction of the superfluid onto the normal fluid. 

The aim of this work is to gain insight into self-consistent turbulent two-fluids hydrody- 
namics using a shell model ij]. In classical fluid dynamics 15| shell models are idealized 
truncated models of the Navier-Stokes equation which neglect any geometry of the flow but 



describe the Richardson cascade, its properties (e.g. the Kolmogorov energy spectrum) over 
a wide range of scales and times. In the less understood context of helium II, shell models 
allow us to explore the interaction of normal fluid and superfiuid in the inertial and dis- 
sipative ranges in a relatively simple way. What happens at very large length scales [l6| 
where finite-size effects may affect the dynamics of the energy- containing eddies is beyond 
the scope of this work. 

II. MODEL 

A. Classical GOY model 

The classical Navier-Stokes equation is 

p(^ + u-Vu) = -Vp + /iV2u, (1) 

with the incompressibility condition V ■ u = 0, where u is the velocity, p the pressure. 



p the density and fi the viscosity. Starting from the works of Gledzer 17| and Yamada & 



Okhitani [l8| , a variety of shell models have been developed [l9[ to study turbulent solutions 
of Eq. [H The GOY model, named after the pioneering authors, is a finite-dimensional 
dynamical system which describes the dynamics of idealised (complex) Fourier components 
of the velocity field, Um corresponding to wavenumber km] the index m {m = 1,- ■ ■ ,M) 
is called the shell index. The wavenumbers are geometrically distributed in k-space, km = 
koX"^, where ko is a reference wavenumber and A > 1. The equations of motion are 

( di ^ ^^^ I ^"^ ^ ^rn[u] + fSm,m', (2) 

for m = 1, ■ ■ ■ M, where u = fi/p is the kinematic viscosity, 6m,m' is Kronecker's delta, and 
/ is the amplitude of external forcing applied on a particular shell m = m'. The inertial 
term Gmfn] is quadratically nonlinear and local in k-space, coupling Um with its nearest 
neighboring shells; more precisely we have 

Gm[u] = «(c^^Um+l^m+2 + c'i^^ Um-lUm+1 + C^^J Um~lUm~2) , (3) 

where u denotes the complex conjugate of u, and 

cW = akm, C^m^ = bkm-l, C^m^ = c/c^-2, (4) 



The boundary conditions are Um = for m < 0, m > M and 

„(2) _ „(3) _ (3) _ (1) _ (1) _ (2) _ ^ /^x 

The coefficients a = 1, b = —1/2, c= —1/2 and A = 2 are chosen so that in the steady 
{d/dt = 0), unforced (/ = 0), inviscid {u = 0) case the nonhnear interaction conserves the 
two quadratic invariants of the 3- dimensional Euler equation, energy and hehcity, 

E=^j^u-udV, H=^jivndV, (6) 

(where a; = V x u and V is volume) which are identified as 

M M 

-E = - ^ \Um\^ = ^ Emkm, (7) 

m=l m=l 

where 

is the spectral energy associated with shell m, and 

1 ^ 

H = -Y,i-irx-w\'. (9) 



m=l 



B. Two-fluids equations 

The two-fluids equations of helium II are [1|] 



O S S 

p'(-^ + u'-Vu') = -^Vp + p'SVT-F (10) 



„n 



p"(^ + u" ■ Vu") = -^V|> - p'SVT + /iV'u" + F (11) 

with V ■ u'^ = V • u" = 0, where S is the specific entropy. Hereafter we denote by i/" = fi/p^ 
the kinematic viscosity of helium II. If F = and V x u" = cj'^ = 0, Eqs. [10] and [TT] 
reduce to Landau's two-fluids equations, which describe the well-known mechanical and 
thermal behaviour of helium II in the absence of vortex lines. In the presence of vortices, 
the mutual friction F which couples the two fluids depends on the vortex line density L 



(vortex length per unit volume), or, more precisely, on some suitably defined coarse-grained 
vorticity field a;*. Unfortunately the form of F is uncertain. The form given by Hall & 
Vinen 20| accounts for experiments in rotating cylinders and Taylor-Couette fiow 2l|, in 
which vortices are either straight or curved, but still polarised in the same direction. In the 
case of turbulent fiow, vortex filaments may be random or partially polarised, so the relation 
between the coarse-grained vorticity cj* and the vortex line density L is not clear. Hereafter 



for the sake of simplicity, we assume the Gorter-Mellink 22| form 

F = ^^nL{u^ - u"), (12) 

where i? is a known temperature dependent mutual friction coefficient. This form for the 
mutual friction is motivated by dimensional arguments and is consistent with the method 
employed to measure vortex line density in experiments. It has been widely used in the 



quantum turbulence literature 
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C. Two-fluids GOY model 

The natural generalization of Eq. |2] to two fiuids is 

^< = Gm[u'] - F„ + rSrn^m': (13) 

where the nonlinear term for each fiuid is as in Eq. [3l and the mutual friction is 

Fm = afi;L«, -u^). (15) 



The temperature dependence of the friction coefficient a = Bp^/{2p) is well-known |2j 
For consistency, the vortex line density L is not an arbitrary parameter, but is identified as 
L = Q}!"^ I K where the superfluid enstrophy Q is 

M 

Q = J2 2^^l<l'- (16) 

m=l 

Clearly this model, which describes superfiuid vorticity as a continuum, is meaningful 
only for length scales bigger than the average intervortex spacing i ^ L~^/^, that is to say 
for wavenumbers km < ki = 1/i = L}^"^. 



In a recent paper, Roche et al. [23| solved equations similar to our Eqs. [TO] and [TT] in the 
presence of continuous forcing applied to the fluid with the greater density, thus performing 
the first Direct Numerical Simulations (DNS) of two-fluids hydrodynamics. Unlike our 
two-fluids shell model, their two-fluids DNS contained an artificial superfluid viscous force 
jjSpS-^2^s ^^ |.]-^g right hand side of Eq. (TDl The introduction of this unphysical term was 
motivated by numerical analysis. They set the ratio u^/u'^ = 4 to minimize the effect of 
the artificial viscosity while preserving the stability of the calculation. We stress that our 
two-fluids shell model (see Eq. [T3|) does not need such an artificial term. Another related 



model is the large-eddy calculation of Merahi et al. |28| 



III. RESULTS 

A. Forced Turbulence 

Unless stated otherwise, the results which we present are obtained by numerically in- 
tegrating in time Eqs. [13] and [14] for M = 18, ko = 2~^. The same forcing (typically 
js _ jn _ (^1 -|- ^) X 5 X 10"'^) is applied to shell m = 4 for both superfluid and nor- 
mal fluid. A typical initial condition is ul^ = (1 -|- i)kmexp (— A;^/2), (j = n, s), for which 
-E'4 = /CmGxp (— /c^). The time stepping combines the Crank-Nicolson method for the dif- 
fusion term and the Adams-Bashforth method for the other terms; the typical time step is 
At = 5 X 10^^. We stop the time integration after the spectrum saturates, making sure that 
the elapsed time is of the order of 10 large eddy turnover times (about 100 turnover times 
of the forcing wavenumber), which is typically of the order of 500 s. We find that, once 
saturation is achieved, time averaged spectra are the same as ensemble averaged spectra 
(obtained using randomly phase shifted reahsations) . 

To study the temperature dependence of the results, we consider three temperatures 
T = 2.157 K, 1.96 K and 1.44 K which hereafter we refer to as high, medium and low 



temperatures 23|]. The values of p^, p", p, z/ and a in TABLE 1. corresponding to these 
temperatures are The resulting normal fluid fractions, p"'/p'^, are approximately 10, 1 and 
0.1 for high, medium and low temperatures respectively. 

Fig. [1] (top) shows superfluid and normal fluid spectra at high temperature. It is appar- 
ent that superfluid and normal fluid velocity lock onto each other by the mutual friction 
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r(K) 


P' (g/cm3) 


pn (g/^^3) 


P (g/cm3) 


u (cm /s) 


a 


2.157 


0.01510 


0.12939 


0.14449 


1.526 X 10"^ 


1.045 


1.96 


0.07335 


0.07221 


0.14556 


9.694 X 10-5 


0.245 


1.44 


0.13251 


0.01264 


0.14515 


9.538 X 10-5 


0.059 



TABLE I. Helium parameters used in our calculations. 



over many length scales as envisaged by Vinen and Niemela |7| and Barenghi et al. 29 |. 
Furthermore a closer inspection of the complex Um throughout the inertial subrange reveals 



that 3? iul 



3fJ iu'L) and ^ [u 



m) 



'^{ut 



to at least 4dp. Both spectra are consistent 
jy the solid line) over a wide inertial range 



with the Kolmogorov scaling k~^^^ (denoted 
ko <^ k <^ ki, as observed in experiments lOj and DNS of two-fluids hydrodynamics 23 |. 
The vertical dotted line denotes the wavenumber k£ which corresponds to the intervortex 
spacing. To make the Kolmogorov scaling more evident, we plot the compensated spectra 
/c^/^ Ek (see Fig. [1] (top, inset)). It is apparent that the superfiuid spectrum extends to bigger 
wavenumbers than the normal fluid's: this is because there are no viscous forces acting on 
the superfiuid; nevertheless, the superfiuid spectrum decays at large k because superfiuid 
motion is damped by mutual friction. We confirmed that the development of a k~^^^ spec- 
trum is dependent neither on our choice of initial condition nor on the fact that we force 
both fluids. 

It is instructive to consider the scale-by-scale energy budget per unit mass. Using Eqs. 
and [H] and the fact that 



we obtain 



dE„ 



dt 



dt 



du 



dt 



m — 



Ur 



~dt 






"inj-- 



rjis 



^m + ^inj^m,4y 



(17) 
(18) 



where 



r^ = -S(G„[n>g (19) 

are the energy transfer rates arising from the triadic interactions between Fourier modes 



within each fluid, 

Dl = -2Z/" klEi:, (20) 

is the rate of viscous dissipation in the normal fluid, 

M;^ = ^5? (i^rnO (21) 

and 

M4 = -3? (F^O (22) 

are the rates of exchange of kinetic energy between the two fluids due to mutual friction 
(M^ being the flow of energy from the normal fluid to the superfluid and M^ the flow from 
the superfluid to the normal fluid) and 

4n, = rSm,m''ii (<) (23) 

are the rate of influx of energy due to the forcing terms, where 9fJ and 3 denote real and 
imaginary parts respectively. The nonlinear triadic interaction T^ is deflned such that the 
energy flux at shell m* is given by 

K= E ^- (24) 

m< m* 

= K+1 - i-b - l)K, (25) 

where we deflne the correlator 

A^„ = kn^iQ «_iu>^„+i) , (26) 

all other terms cancelling due to the conservation of energy. 

Fig. [1] (bottom) shows the energy balance in the inertial range. It is apparent that the 
total energy flux in each shell is zero (solid black diamond). We flnd that the normal fluid's 
inertial term (hollow red squares) is balanced by the viscous term (blue hollow triangles), 
as in ordinary turbulence. On the contrary, the superfluid's inertial term (solid red squares) 



is balanced by the mutual friction (solid grey circles), in agreement with Roche et al. 23 |. 
The inset of Fig. [1] (bottom) shows the energy balance over the entire k-space, including the 
contribution of the forcing on the shell m = 4. We observe that M^ and M^ are always of 
opposite signs, as expected from their definitions, and, almost without exception, M^ > 
and M^ < 0. Furthermore D'l^ < for all m. 



The basic picture did not change when we reduced the temperature to T = 1.96 K 



(Fig. |2]), and then to T = 1.44 K (Fig. [3]), again in agreement with experiments 10|. There 
are however noticeable differences which develop as the temperature is lowered. At high 
temperature viscous forces provide a sufficient energy sink for both normal fluid and su- 
perfluid. As the temperature drops the normal fluid fraction decreases and the superfluid 
fraction increases. This reduces the ability of the normal fluid viscosity to adequately dis- 
sipate the energy. As T is decreased, both the normal fluid and superfluid spectra extend 
to the right, towards higher wavenumbers, and the wavenumber, ke, corresponding to the 
intervortex spacing, moves progressively to the right. Furthermore at low T in the normal 
fluid the viscous dissipation is no longer balanced by the inertial term, but by the mutual 
friction term (see Fig. [3] bottom). 

In their pure states {p'^ = and p" = respectively) we expect the spectrum of the 
normal fluid to decay exponentially after k^ (see Eq. I7r|) . the Kolmogorov microscale at 
which the inertial term is balanced by the viscous dissipation, and the superfluid to continue 
to cascade like k~^'^; a build-up of energy at high k is thus expected, although the physical 
interpretation of our model is limited to k < kg. We can furthermore expect that as we 
approach these limiting cases the dominant fluid will cause a deviation in the spectrum of 



the other fluid. L'vov et al. |30| derived approximate expressions for each of these cases. 



At low temperatures (p* ^ p""), the normal fluid spectrum, instead of decaying exponen- 
tially for k > kr^, deviates slightly from k~^^^ due to the force exerted on it by the mutual 
friction. This continues until the length scale, denoted k^ (see Eq. |28|) . at which the viscous 
dissipation balances the mutual friction. Beyond this wavenumber the normal fluid spectrum 
decays with a power-law A;~^^/^; not exponentially as in the case of classical turbulence. This 
power-law decay is valid at least until ki, the wavenumber corresponding to the inter-vortex 
spacing, at which point the model for the superfluid is no longer valid. 

On the other hand, at high temperatures (p" ^ p*), the superfluid spectrum beyond 
k^i is affected by the exponentially decaying normal fluid so that it deviates from its k~^^^ 
power-law to a steeper slope k~^. This continues for as long as the mutual friction dominates 
over the superfluid inertial term. At some wavenumber, k+ (see Eq. I29l) . this relationship 
shifts such that the inertial term becomes dominant and the superfluid spectrum regains the 
/;;-5/3 power- law for as long as the model remains valid {k < ki). 

9 



Following LVov 3C|, we define 

K = ^"'^-"\ (27) 

/ s \ 1/2 

1/2/^1 nV4,-i/2 



K = a'" ^-j Q^/V-V^ (28) 

fc+ = a3/2g3/4^-i/2^ (29) 

where v = v"^ and e = —dE/dt is the mean energy dissipation (per unit mass), which, in 
the inertial subrange, is approximately constant. In the context of shell models we define 

We furthermore find, due to the efficient locking of the superfluid and normal fluid velocities 
in the inertial subrange, that e" ~ e*, which we refer to as e in Eqs. f l27|29|) . 

We show the deviation at low temperature in (Fig. H]), in which we show both power-laws 
and the critical wavenumbers at which the spectra change. We calculate these wavenumbers 



in accordance with the approximations derived by L'vov et al. 30|. The deviations at high 
temperature are best realised by considering a system of turbulent superfluid ^He-B (see 
below) . 

As mentioned above, we expect that at a sufficiently low temperature the superfluid 
energy spectrum builds up at high k. This is because the only energy sink for the superfluid 
is the mutual friction which depends on the presence of the normal fluid to dissipate the 
superfluid energy. We show the development of the build-up of energy in Fig. [5l We 
produced Fig. [5]by allowing the fluids to attain fully-developed spectra at high temperature. 
We then reduced the temperature significantly to a finite, fixed temperature well below our 
previous low T and monitored the changes in the superfluid spectrum over time. In order to 
maintain numerical stability, we fixed the values of p*, p" and z/". Thus the only temperature- 
dependent variable was a, which alone dictates the magnitude of the mutual friction felt 
by the superfluid. The spectrum is seen to change as progressively more energy moves to 
higher k. This phenomenon highlights the necessity for an energy sink at high k such as the 
Kelvin wave cascade and phonon emission. 

In classical turbulence the Reynolds number. Re = UD/u, is a measure of the ratio of 
the inertial and viscous terms, where U and D are the large scale velocity and length scale 

10 



respectively. The inertial range exists as long as Re ^ 1. In the same spirit we may define 
a quantum Reynolds number as the ratio of the inertial and general dissipative terms, be 
they due to viscous dissipation or mutual friction. Using this definition we may say that 
the Reynolds number for the superfiuid is temperature-dependent, in agreement with 23| : 
both normal fiuid and superfiuid spectra extend to higher wavenumbers as the temperature 
is decreased. 

We furthermore show the dependence of the vortex line density, L, on the Reynolds 
number. According to Kolmogorov, the ratio of the dissipation scale and the scale of the 
large eddies is 5/D oc i?e~^/^. In superfiuid turbulence the smallest scale is £ oc L~^/^, 
hence we expect L oc Re^'"^. We show that this scaling holds true by considering Re and 
L and allowing the forcing to vary between / = (1 + i) ■ 5 ■ 10^^ and / = (1 + i) ■ 5 ■ 10^ in 



mu 

al. 



tiples of 10 (see Fig. [6]). This finding agrees with the truncated DNS model of Salort et 



m 



B. Decaying turbulence 

To study the decay of turbulence we let / = and start from saturated spectra as an 
initial condition. We observe that during the decay both the superfiuid and normal fiuid 
spectra maintain their initial shape {k~^^^ at low to intermediate values of k, followed by a 
more rapid drop at larger k), as shown in Fig. [7] (top and bottom). The spectra shown in 
these figures are the result of ensemble averaging over 10 realisations. We also found that 
the total turbulent kinetic energy, E{t), where 

E\t) = J2l\uUt)\', {J = n,s), (31) 

m 

and the vortex line density, L(t), decay as E(t) oc t~^ and L(t) oc t~^^'^ as shown in Figs [8] 
and[9l which is in agreement with experiments 32|-|34| and theoretical models 35 1. 



C. Forced turbulence in superfiuid "^He-B 

The other non-radioactive isotope of liquid helium is ^He. ^He atoms are fermions, having 
only a single neutron in their nuclei, and achieve superfiuidity by the Cooper pairing of two 
atoms. Although T^ for ^He is about 2.17 K, for ^He it is much lower, at about 1 mK. 

11 



Experimental studies of ^He have focused on the so-called B-phase, known as ^He-B. The 
most notable physical property of '^He-B is the rapid increase in the normal fluid viscosity 



below Tc, which is proportional to 1/T^ J36l]. In fact the viscosity is so great that for 



all intents and purposes the normal fluid may be considered to be at rest. We may thus 
consider such a system to be that of a turbulent superfluid in the presence of a stationary 
normal fluid. This description is equally applicable to that of high temperature ^He in the 
wayenumber subrange k ^ kj^. This correspondence has already been noted by LVov et 



al. [30|. A second relevant physical property is the change in its quantum of circulation such 
that K = /i/2m3, or about 2/3 that of '^He, where m^ denotes the mass of a '^He atom. 



In a theoretical paper Vinen 26| showed that in superfluid ^He-B the positions of the 
dissipative subrange and the inertial subrange are reversed, such that the dissipation due 
to the mutual friction occurs at low wavenumbers whereas the inertial subrange obeying 
the /c~^/^ power-law is located at high wavenumbers. His prediction for the wavenumber at 
which this changeover takes place is in agreement with that of L'vov et al. J30|, |37|. LVov et 
al. [37| furthermore predicted that the dissipation should follow a k~^ power-law, the same 
power law as that of high temperature ^He in the wavenumber subrange k ':$> krj (30|]. 

In order to model turbulence in '^He-B we reduced our two-fluids shell model back to a 
single-fluid model, replacing viscous dissipation with a dissipation due to mutual friction. 
The shell model equation then becomes 

where the mutual friction is now of the form 



F^ = aK ""^Lu'^. (33) 

The temperature dependence of Eq. [32] is controlled by the parameter a alone. Experience 
has shown that a mutual friction of this form is extremely efficient at dissipating energy 

nnn 

and in order to attain a steady spectrum of the kind described in [26|, |30|, |37[ it is necessary 
both to reduce a to O(10~^) (much below Tc and very close to a pure superfluid state) and 
to increase the forcing to O(IO^) (5 orders of magnitude greater than that used for ^He). 
The numerical results show the spectrum initially decreasing like k~^ which subsequently 
changes to k~^/^. The changeover wavenumber is in reasonably good agreement with fc_|_ 
predicted by L'vov (see Fig. [10]). We used the following parameters (M, A, m') = (20, 2, 4). 
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IV. DISCUSSION 

We have developed a two-fluids shell model based on the GOY shell model of classical 
turbulence. We have shown that, in continually excited turbulence, the mutual friction ef- 
fectively locks together the normal fluid and superfluid over a wide range of wavenumbers, 
forming a k~^'^ Kolmogorov spectrum. By reducing the temperature and hence the normal 
fluid fraction, we have shown that this range's cut-offs are temperature-dependent for both 
the normal fluid and superfluid. In other words, lowering the temperature leads to longer 



inertia, 
ing 



subranges. This result for the energy spectrum agrees with the experiment of Tabel- 

. The relation which we find between vortex line 



10| and the DNS of Roche et al. 



density and Reynolds number agrees with Salort et al. 3l| . We have realised the deviations 
from the /c"^'"^ spectrum due to the mutual friction predicted by L'vov 30| for ^He at low 
temperature and by Vinen 26| and L'vov |30l . |37| for ^/fe-B. In the zero-temperature limit 
our model develops a build-up of energy at high wavenumbers, highlighting the necessity for 
an energy sink at high fc, such as the Kelvin wave cascade and phonon emission 



. Our 



model shares the common advantage of all shell models: being flexible and computationally 
inexpensive, it is relatively easy to investigate the decay of turbulence, which would be more 
difficult with DNS. Our results for the decay of turbulence are in agreement with experi- 



ments 



32 



-[3J] and theoretical models 35|. We anticipate that our two- fluids model can be 



used to study other aspects of quantum turbulence. 
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FIG. 1. (Colour online). Top: Main plot: Log-log plot of superfluid (blue diamond) and normal 
fluid (red circle) energy spectra E^ (cm^ s~^) vs wavenumber k (cm~^) at T = 2.157 K. The 
solid line denotes the k~^'^ Kolmogorov spectrum. The vertical dotted line marks k^. Inset: 
Compensated spectra. As above, but k^'^E^ vs k. Bottom: Main plot: Log-lin plots of time- 
averaged scale-by-scale energy budget dE^/dt (cm^ s~^) vs wavenumber k (cm"-*^) for the two-fluid 
model at high temperature (2.157K). We show the fluxes for the normal fluid (hollow shapes) and 
superfluid (solid shapes). The main plot shows the balance of fluxes in the normal fluid between 
the inertial term T^ (red squares) and the viscous term D"^ (blue triangles), and in the superfluid 
between the inertial term T^ (red squares) and the mutual friction term M^ (grey circles). We 
also show the total flux dEm/dt (black diamonds). Inset: In the inset we show the fluxes over the 
entire wavenumber range, including the flux due to the external forcing e^^- (light blue line) which 
is nonzero only at shell n = 4. The vertical dotted line denotes the time-averaged wavenumber, k^, 
corresponding to the intervortex spacing. 
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FIG. 2. (Colour online). Top: Main plot: Log-log plot of superfluid (blue diamond) and normal 

fluid (red circle) energy spectra E^ (cm^ s~^) vs wavenumber k (cm~^) as in [1] (top) but at 

T = 1.96 K. The solid line denotes the k~^'^ Kolmogorov spectrum. The vertical dotted line 

marks kg. Inset: Compensated spectra. As above, but k^'^Ek vs k. Bottom: Main plot: Log-lin 

plots of time-averaged scale-by-scale energy budget dEj^/dt (cm^ s~^) vs wavenumber k (cm~^) for 

the two-fluid model as in[T] (bottom) but at medium temperature (1.96K). We show the fluxes for 

the normal fluid (hollow shapes) and superfluid (solid shapes). The main plot shows the balance 

of fluxes in the normal fluid between the inertial term T^ (red squares) and the viscous term D^ 

(blue triangles), and in the superfluid between the inertial term T^ (red squares) and the mutual 

friction term M^ (grey circles). We also show the total flux dEm/dt (black diamonds). Inset: 

In the inset we show the fluxes over the entire wavenumber range, including the flux due to the 

external forcing e^^- (light blue line) which is nonzero only at shell n = 4. The vertical dotted line 

denotes the time-averaged wavenumber, ki, corresponding to the intervortex spacing. 
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FIG. 3. (Colour online). Top: Main plot: Log- log plot of superfluid (blue diamond) and normal 
fluid (red circle) energy spectra E^ (cm^ s~^) vs wavenumber k (cm~^) as in [1] (top) but at 
T = 1.44 K. The solid line denotes the k~^'^ Kolmogorov spectrum. The vertical dotted line 
marks ki. Inset: Compensated spectra. As above, but k^'^Ek vs k. Bottom: Main plot: Log-lin 
plots of time-averaged scale-by-scale energy budget dEj^/dt (cm^ s~^) vs wavenumber k (cm~-^) 
for the two-fluid model as in[T] (bottom) but at low temperature (1.44K). We show the fluxes for 
the normal fluid (hollow shapes) and superfluid (solid shapes). The main plot shows the balance 
of fluxes in the normal fluid between the mutual friction term M"^ (red squares) and the viscous 
term DJ^ (blue triangles), and in the superfluid between the inertial term T^ (red squares) and 
the mutual friction term M^ (grey circles). We also show the total flux dEm/dt (black diamonds). 
Inset: In the inset we show the fluxes over the entire wavenumber range, including the flux due to 
the external forcing e^^- (light blue line) which is nonzero only at shell n = 4. The vertical dotted 
line denotes the time-averaged wavenumber, kp, corresponding to the intervortex spacing. 
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FIG. 4. (Colour online). Log- log plot showing dissipation of turbulent energy in normal fluid at 
low temperature. Coupled normal fluid (red circles) is seen to deviate from uncoupled normal fluid 
(green triangles) as a result of mutual friction with superfluid (blue diamonds). Coupled normal 
fluid follows k~^'^ for k < A;^, deviates slightly for k^ < k < k^, and follows a A;~^^'^ power-law for 
k > k^. The vertical dotted line is k^ and the short- and long-dashed lines are respectively kr^ and 
ki,. For clarity krj < k-^ < ki. 
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FIG. 5. (Colour online). Log- log plot of development of build-up of energy in superfluid spectrum. 
Spectra (bottom to top) at time t = (red), 0.025 (green), 0.05 (dark blue), 0.5 (pink), 5 (light 
blue), 50 (black), 500 (orange), 5000 (grey) and 10000s (red) after lowering the temperature, ki 
(not shown) moves from k ~ 10^ to /c ~ 2 • 10^ during this period. The thick black line denotes 
the k~^i^ Kolmogorov spectrum. 
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FIG. 6. (Colour online). Log-log plot of vortex line density, L, vs Reynolds Number, iZe, at all three 
temperatures: high (red squares), medium (grey circles) and low (blue triangles). The points from 
left to right correspond to increasing forcing, the leftmost point using forcing / = (1 + i) • 5 • 10"^ 
and the rightmost point / = (1 + i) • 5 • 10^. The Reynolds Number is calculated at the first shell 
(n = 1) giving D = 2^. The solid black line is L = Re^'"^ . 
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FIG. 7. (Colour online). Decay of energy spectra for normal fluid (top) and superfluid (bottom) 
over 10 realisations. Spectra top to bottom: After 500, 1000 , 2500 and 5000s. Also shown is the 
^-5/3 spectrum (solid black line). 
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FIG. 8. (Colour online). Decay of total energy over period of 5000s ensemble averaged over 10 
realisations. Top to bottom: High (red line), medium (grey line) and low (blue line) temperatures. 
Shifted to show power law. Also shown is the i~^ spectrum (solid black line). 
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FIG. 9. (Colour online). Decay of vortex line density over period of 5000s ensemble averaged 
over 10 realisations. Top to bottom: Low (blue line), medium (grey line) and high (red line) 
temperatures. Data is not shifted. Also shown is the i~^'^ spectrum (solid black line). 
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FIG. 10. (Colour online). Log-log plot of fully developed, steady power spectrum for superfluid 
(red circles) in the presence of a stationary normal fluid in '^He-B. We show the two power-laws, 
k~^ and k~^'^ . The dashed line is kj^, the wavenumber at which the inertial term becomes of the 
same order as the mutual friction term, and the dotted line is /c£, the intervortex spacing. 
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